%y是当前所在列，x是将填入行数-1,J是病种数，N是病人数，n是剩余病人个数，j是剩余病种数
function[A]=findPIE2(A,J,j,n,y,x)
if y<J
    for i=0:n
    x0=x+1;x=x+nchoosek(j-2+n-i,j-2);%填入数据的首行与尾行
    A(x0:x,y)=i;
    A(x0:x,J)=A(x0:x,J)-i;%填入数据并修改最后一列数据
    [A]=findPIE2(A,J,j-1,A(x0,J),y+1,x0-1);%对下一列应用该程序进行修改，套娃      
    end
end
end